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Abstract 

We consider a classical toy model of a massive scalar field in 1 + 1 dimensions with a constant ex- 
ponential expansion rate of space. The nonlinear theory under consideration supports approximate 
oscillon solutions, but they eventually decay due to their coupling to the expanding background. 
Although all the parameters of the theory and the oscillon energies are of order one in units of the 
scalar field mass m, the oscillon lifetime is exponentially large in these natural units. For typical 
values of the parameters, we see oscillon lifetimes scaling approximately as r oc e~x.p(kE/m)/m 
CN . where E is the oscillon energy and the constant k is on the order of 5 to 15 for expansion rates 

between H = 0.02m and H = 0.01m. 
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INTRODUCTION 

Central to many unsolved problems in particle physics and cosmology, ranging from 
the flatness of the slow-roll inflaton potential to the hierarchy of scales in grand unified 
theories, is the need to understand the origin of "unnaturally" large or small dimensionless 
parameters. In this Letter we demonstrate how oscillons — localized, oscillatory solutions to 
nonlinear field equations that either are stable or decay only after many cycles — can provide 
an explicit example in which such unnatural behavior emerges dynamically. We show that 
in the presence of an expanding background metric, oscillons in a simple toy model have 
lifetimes that are exponentially large compared to the natural scales of the system, even 
though all the inputs to the theory are of order one in these units. Although the situation 
is reminiscent of tunneling behavior in quantum systems, this decay is entirely classical. 

A wide variety of nonlinear field theories have been found to support oscillon solutions 
(also known as breathers) ll|, |2J |3|, |4J, la, laM, la |9|, |lC| |ll(. In some special cases, such as 
the sine- Gordon breather |1[ and Q-ball |2|, conserved charges guarantee the existence of 
exact, periodic solutions. In <p A theory in three dimensions, in contrast, y,|5(, oscillons are 
observed to decay suddenly after lifetimes of order 10 3 to 10 4 in natural units. In other 



three-dimensional models 11], no sign of decay is visible in current numerical simulations. 
In a broad class of one-dimensional models, perturbative analyses [l| and numerical sim- 
ulations J3, |9| suggest that the oscillon lifetime is infinite, but analytic arguments |l2j point 
to the existence of non-perturbative, exponentially suppressed decay modes. Such models 
then have a curious property: in a simple, classical theory in which all parameters are of 
order one in natural units, classical dynamics generate a quantity — the oscillon lifetime - 
that is exponentially large compared to the natural scales of the system. 



MODEL 

In this Letter, we investigate explicitly a simple example of such a system. We consider a 
toy model consisting of a classical real scalar field <f>(x, t) of mass m in one space dimension, 
governed by the Lagrangian density 
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In this classical theory, the dynamics are unchanged when the Lagrangian density is multi- 
plied by a constant (equivalent to rescaling H). We have used this freedom, combined with a 
rescaling of the field 0, to fix the magnitude of the coefficient of the </> 4 term while maintain- 
ing the conventional form of the free term. The arbitrary choice of the ratio of magnitudes 
of the 4 and 6 terms is chosen for convenience and is not essential to our results. In the 
natural units of the mass m, then, all the coefficients of the potential are of order 1. The 
only important choices are the signs: the sign of the <fi 2 term gives the <fi field a conventional 
(non-tachyonic) mass; the sign of the 4 term is necessary for oscillons to exist; and the 
sign of the 4> 6 term is necessary for the field to be stable against runaway growth. We note 
that we have observed the same behavior in 4 theory with the standard symmetry-breaking 
potential. Here we have chosen the 6 model just to avoid any distractions that might be 
caused by the existence of static "kink" solutions in the 4 theory; our theory contains no 
static solutions. 



Like the other one- dimensional models discussed above, in a static background the model 
of Eq. (£Q) supports oscillons that appear to live indefinitely in numerical simulations. They 
are spatially localized, with size of order 1/m and fundamental frequencies of oscillation 
comparable to but always below m. We would like to include an additional element: coupling 
to an expanding background metric, inspired by the consideration of oscillons in the early 
universe (work in pro gress 13l ; see also studies of nontopological solitons and oscillons in 
hybrid inflaton [1J, [l5( and axion-oscillons |l6[ , and for a broader review of lattice field theory 
simulations in the early universe see |l7|). In comoving coordinates, then, the Lagrangian 
takes the form 
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leading to the equation of motion 
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where the physical coordinate is related to the comoving coordinate x by the scale factor 
a(t) and we hold the Hubble constant H = a(t)/a(t) fixed, giving an exponential expansion 
rate. Here <fi is the derivative of <fi with respect to time and <$' is the derivative of <fi with 
respect to the comoving coordinate x. This expansion destabilizes oscillons whose width is 
of order of the horizon size 1/H, but in numerical simulations we observe that oscillons of 
smaller size remain stable, maintaining a fixed size in physical units. 

A technical advantage of this model is we can study oscillon behavior efficiently for 
extremely long times, because any regions at distances significantly greater than the horizon 
length 1/H from an oscillon cannot influence its evolution. Thus we do not need to expand 
the box proportionally to the runtime [ll| or introduce absorbing boundaries |l8( in order 
to prevent unwanted reflections from disturbing the oscillon solution under study. 



NUMERICAL SIMULATION 



Starting from thermal initial conditions with T > m, we see that oscillons emerge co- 
piously as the expansion cools our universe. These initial conditions are generated using 
a canonical ensemble of quantized modes of the free scalar field, and thus contain no fine- 
tuning. 1 To simplify the classical analysis, we neglect zero-point fluctuations. Oscillons have 
been similarly found to emerge spontaneously during phase transitions |15|, |l9( . 



1 Here we must do a quantum calculation to avoid the Jeans paradox. As a result, h appears explicitly and we 
can no longer scale out the overall normalization of the coupling constants through our choice of units, as we 
have in the classical calculation above. The equal-time commutation relation [(f)(x, t), <j)(y, t)\ = ihS(x — y) 
implies that 4> has units of \fh. Therefore the Lagrangian density takes the form 
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where g is a constant with units of 1/H. In the classical theory, we may choose our units so this quantity is 
equal to 1. In the thermal calculation, however, we encounter the quantity hu/T in the Boltzmann factor 
for modes of energy lo. Choosing units so that ft = 1 in this expression forces us to include g explicitly in 



We use a straightforward numerical simulation in which we discretize x at the level of 
the Lagrangian in Eq. (J2J), working in natural units where m = 1. For the space derivatives 
we use ordinary first-order differences, 

0W) = <W|zM) (5) 

where cj) n refers to the value of at lattice point n. We work on a regular lattice with spacing 
Ax = 0.01, and impose periodic boundary conditions. Varying this Lagrangian yields lattice 
equations of motion with second-order space derivatives, 

7 W ,^); /,x 0n+l(*)+0n-l(*) -20 w (t) ..3 5 

<M*) + ^y<W) = a(t) 2 Ax 2 <Pn(t) + <p n (t) - (p n (t) . (6) 

Finally, we use second-order differences 

7 /.x __ <Pn{t + At) + w (t - At) - 20„(t) ^ . __ w (t + At)-0 w (t-At) 

with At = 0.005 to compute the field at t + At based on the values at t and t — At. We 
have verified that our results are not sensitive to the particular choice of lattice spacing and 
time step. (By the Courant condition we must choose At < Ax or the algorithm will be 
unstable.) 

In the absence of expansion, our system conserves energy exactly in the limit of infinites- 
imal time step At (regardless of the spatial lattice size). In an expanding background, the 
configuration loses energy at a rate given by the pressure times the rate of change in volume 

— — — — p(x,t)a(t)dx (8) 



where the pressure density is 
P(x,t) = - 
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We use this result to check the accuracy of our calculation, and find that it is maintained 
to better than 1 part in 10 4 , as we'd expect for a second-order method. As the universe 
expands, our lattice expands with it, but the oscillons do not. Therefore, whenever our lattice 
has expanded by a factor of 2, we refine the lattice, doubling the total number of lattice 
points and bringing the lattice spacing back to its original value in physical units. We assign 
values to the field for the new intermediate lattice points by polynomial interpolation. This 



the classical Lagrangian density. An oscillon solution with g = 1 will continue to be a solution for general 
g, provided that we scale up its amplitude by a factor of 1/y/Kg. For our classical treatment of the time 
evolution to be valid, the oscillon must have an energy that is large compared to the energy of a typical 
quantum fluctuation fiwo, where luq ~ m/h is its fundamental frequency of oscillation. Since the oscillon 
energy goes like m/(Hg), this requirement can be implemented by choosing kg to be small — the classical 
approximation is valid at small coupling. We have verified that our results are not affected by the value 
of this parameter; as long as T> m/(Hg), oscillons are generated copiously. They then evolve according 
to the same classical dynamics with the rescaled field. 



numerical algorithm is highly stable, maintaining precision without any sign of degradation 
even after extremely long runs. 

We begin our simulation using values of <f> n (t = 0) and <fi n (t = 0) drawn from a random 
thermal configuration with T = 1 (setting the quantum scaling parameter g to one for 
simplicity). As the expanding universe redshifts away the ordinary perturbative fluctuations, 
it becomes simple to pick out the oscillon peaks. At this point, we excise a particular 
oscillon, keeping only a window of size 100 + 2/H in physical units around it. (We identify 
the starting and ending points of the oscillon profile as the places where the derivative of 
the energy density is 0.001 in natural units.) Now, each time we insert new lattice points, 
we also truncate the lattice back down to this size, so that the total computational cost of 
the run scales only linearly with time. Because of the expansion, any noise introduced by 
this truncation can never affect the oscillon dynamics. We have verified that changing the 
box size used for this truncation does not affect the oscillon lifetimes we observe. 



DECAY ANALYSIS 

Although the oscillon profiles we obtain by this method vary, they all exhibit common 
behavior. We use a small-amplitude asymptotic analysis in a static background, following 
|il,[20j, to gain some insight into the general properties of oscillon solutions. At large \x\, the 
magnitude of 4> is small and we can ignore the nonlinear terms in the equations of motion, 
so our solution must be of the form exp(— m|x|e) smut, where u = m\J\ — e 2 . (We have 
absorbed an arbitrary phase by our choice of the zero of t.) Therefore it makes sense to work 
in terms of rescaled variables x = mxe an d t = mt^Jl — e 2 , giving the equation of motion 



(0) TT + <f> = e 2 \U>)„ + (0) xx j + ^ - <f . (10) 

Next we expand <fi in powers of e using these variables, 

<J){X, r) = e<Mx, r ) + e2 Mx, r) + e 3 3 (x, r) + . . . (11) 

and consider the equations of motion order by order in e. At 0(e), we have simply 

(0i) rT + 0! = O (12) 

and thus 4>\ — /(x) srn T, where we have again absorbed a phase into the definition of r, 
and the profile f(x) remains to be determined. At 0(e 3 ), we have 

(fo) TT + h = (0O TT + (0x) xx + 0? . (13) 

The key point is that for all Xi the function of r on the left-hand side of this equation is 
orthogonal to sinr, so the right-hand side must be as well. As a result, setting the Fourier 
coefficient of sin r on the right-hand side equal to zero gives 

0-/W + !j/(x) 3 = O (14) 

whose solution, consistent with the boundary conditions at infinity, is /(%) = \/| sech^. 
Thus we have a family of approximate solutions 
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parametrized by the small amplitude e. 

Although this analysis is modified by higher orders in e and the inclusion of the expanding 
background, we find numerically that its general features remain. In particular, we can have 
oscillons of large width as long as they have correspondingly small amplitudes; the oscillon's 
total energy, which is proportional to its width times its amplitude squared, scales linearly 
with its amplitude and inversely with its width. Each oscillon decays by a gradual process 
of expansion, through which its width increases and its amplitude and energy decrease. As 
the width approaches the horizon size 1/H, this decay becomes very rapid and the oscillon 
decays suddenly. This process is shown in Fig. ^ We define the decay of the oscillon as the 
time when the derivative of its energy is everywhere below 0.001 in natural units. Since the 
decay process is so sudden, choosing any other reasonable criterion would make a negligible 
change in the lifetime we measure. 

Oscillon Energy vs Time 
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FIG. 1: The final decay of an oscillon as it approaches the point of instability. Although the decay 
is very sudden, the approach to this decay at earlier times is extremely gradual. Here the energy 
is measured in units of m and the time is measured in units of 1/m. 



The dramatic behavior shown in Fig. ^ is visible only at the end of the oscillon's life, 
however. For slightly higher energies, this curve becomes extremely close to horizontal; 
oscillons with only marginally higher energy lie far to the left. As a result, the oscillon 
lifetime scales exponentially with its energy, as shown in Fig. |21 Here we have used the 
energy to select oscillons that we predict will decay in a reasonable amount of time. However, 
we have also verified that more energetic oscillons live for the exponentially long times that 
this analysis would suggest. In particular, we have followed an individual oscillon for times 
exceeding 2 x 10 7 in natural units. This analysis also verifies that the decay is not the result 
of accumulated numerical error. 

To track the flow of energy during the oscillon's decay, we consider a fiducial box with a 



Oscillon lifetime vs Energy 
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FIG. 2: Logarithmic plot of oscillon lifetime as function of the initial 
oscillon energy for various values of the Hubble constant H. The vertical 
axis is the natural logarithm of the lifetime of the oscillon and the 
horizontal axis is the oscillon energy, both in units where m = 1. We 
see that the oscillon lifetime scales like exp(kE/rn)/m, where E is the 
oscillon energy. The rate of exponential growth k is given in as a function 
of H in the table at right. The values of k are roughly linear in H, with 
a slope of approximately — 994/m. 

radius of one Hubble length l/H. The rate of change in the energy within the box is then 
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where the first term represents the flow of of any outgoing waves through the boundary and 
the second represents the energy lost to gravitational expansion. For our oscillon configu- 
rations, both terms oscillate. On average, however, we find the surprising result that the 
first term is both positive and much smaller in magnitude than the second term, meaning 
that the oscillon decays due to energy lost through the expansion rather than through the 
emission of scalar field waves. 

Finally, we note that although the expansion of the universe explicitly introduces an 
exponentially large scale — the size of the universe — it is irrelevant to the oscillon lifetime. 
The oscillon only sees its local region of space, which is expanding at a rate of order one 
in natural units, with greater expansion of the universe corresponding to shorter oscillon 
lifetimes. 



CONCLUSIONS 

We have found that oscillon lifetimes in an expanding background scale exponentially with 
the oscillon's energy in a simple toy model. The oscillon configurations form generically from 
an uncorrelated thermal background. As a result, although all the parameters of our theory 
are of order one in natural units, the dynamics of the theory generates an exponentially large 
time scale through the oscillon's decay. Although our analysis is limited to one dimension, 
we have also conducted preliminary experiments showing similar behavior when the three- 
dimensional oscillons of Refs. p, lll| are coupled to an expanding background in an ansatz 
with spherical symmetry. Work is in progress [13[ to explore the possible consequences of 
these ideas in the early universe. 
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